CaFe2As2: a springboard to investigating Fe-pnictide superconductivity 
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We present detailed electronic structure calculations for CaFe2As2. We investigate in particular 
the 'collapsed' tetragonal and orthorhombic regions of the temperature-pressure phase diagram 
and find properties that distinguish CaFe2As2 from other Fe-pnictide compounds. In contrast to 
the tetragonal phase of other Fe-pnictides the electronic structure in the 'collapsed' tetragonal 
phase of CaFe2As2 is found to be strongly 3D. We discuss the influence of these properties on 
the formation of superconductivity and in particular we find evidence that both magnetic and 
lattice interactions may be important to the formation of superconductivity. We also flnd that the 
Local Spin Density Approximation is able to accurately predict the ordering moment in the low 
temperature orthorhombic phase. 

PACS numbers: 74.70.Dd, 71. 15. Mb, 74.25Jb 
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I. INTRODUCTION 

The recent discovery of superconductivity in the doped 
iron pnictide compounds and subsequent improvement in 
Tc has generated significant interest in uncovering the 
mechanisms responsible for this novel superconductiv- 
ity. In the 1,2,2 class of Fe-pnictide compounds including 
(Ca,Ba,Sr,Eu)Fe2As2, superconductivity has been shown 
under pressure tuning for (Ba,Sr,Eu)Fe2 Asjii^ with tran- 
sition temperatures as high as 32K. There have been 
several reports of the presence of superconductivity in 
the phase diagram of CaFe2As2 as it is pressure tuned 
from the orthorhombic phase into the collapsed tetrag- 
onal phase at low temperatures^!^. The determination 
of the presence of superconductivity has recently been 
shown to be more complex by measurements made using 
a Helium pressure mediuni^. It has also been suggested 
that superconductivity results from the presence of a 
mixed phase intermediate between the collapsed tetrago- 
nal and orthorhombic phases^. Therefore, CaFe2As2 may 
be a material that is near the border for the formation 
of superconductivity and as a result its electronic and 
magnetic structure have the potential to provide criti- 
cal information regarding the superconductivity of these 
compounds. 

The Temperature-Pressure phase diagram for 
CaFe2As2 is unusual amongst Fe-pnictide compounds 
and is shown schematically in Fig. [T] The presence of 
the 'collapsed' tetragonal phase at low temperatures 
with the application of modest pressures distinguishes 
this compound from other members of the 1,2,2 family^. 
The low temperature orthorhombic antiferromagnetic 
phase that occurs at low pressures is common to the 
1,2,2 class and is thought to result from the formation of 
spin density wave (SDW) itinerant magnetism^. At high 
temperatures CaFe2As2 exists in a tetragonal structure 
and is non-magnetio^. 

In the calculations reported here we seek to consider 
the electronic structure of CaFe2 As2throughout the three 
key parts of the phase diagram. We investigate the mech- 
anisms for the novel phase transitions and superconduc- 
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FIG. 1: (Color online) Phase diagram of CaFe2As2 showing 
the Local Spin Density Approximation Fermi surface calcu- 
lated in this study for each region. The dramatic changes in 
the Fermi surface as we introduce new order by structural or 
magnetic phase changes illustrates coupling to the electronic 
degrees of freedom. This is indicative of the need to con- 
sider both magnetic and lattice interactions with the crystal 
structure when discussing superconductivity in this class of 
compounds. 



tivity of this compound. We also present predictions for 
future quantum oscillation studies that may experimen- 
tally verify the electronic structure. 



II. STRUCTURE AND METHOD 

The electronic structures were determined by the Full- 
Potential LAPW method implemented in WIEN2Ki2,. In 
all cases RJ^^Kmax — 8 was used. We have used the 
Local Spin Density Approximation (LSDA) for the cor- 
relation functional except where otherwise stated. We 
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have used the experimental lattice parameters reported 
for each part of the phase diagram. Further details will 
be discussed as is relevant to each calculation. Fermi 
surfaces were visualized with XCrysDenl^ . 

In the following sections we consider each part of the 
phase diagram in turn. First, we present the results for 
the ambient pressure tetragonal phase. These results 
place CaFe2As2 in a background from which it may be 
referenced to other Fe-pnictide systems. We then con- 
sider the collapsed tetragonal phase at P=0.63 GPa and 
contrast this electronic structure with the high tempera- 
ture tetragonal phase. Finally, we consider the antiferro- 
magnetic orthorhombic state and investigate the predic- 
tion of the magnetic moment in this phase. With each 
of these elements of the phase diagram in place, we are 
then in a position to qualitatively discuss their interplay 
in determining the magnetic order and possible mecha- 
nisms for superconductivity in the system. 



III. THE HIGH TEMPERATURE 
TETRAGONAL PHASE 

The crystal structure of the high temperature tetrago- 
nal phase of CaFe2As2 is in the lA/mmm tetragonal 
space group with lattice parameters a — 3.912A and 
b = 11.667AA2,. Calculations of the electronic structure of 
CaFe2As2 in this phase have previously been presented 
for comparison with X-ray photoelectron spectra—. 
Here, we present details, in particular the form of the 
Fermi surface, that are important as a point of com- 
parison between CaFe2As2 and other Fe-pnictide sys- 
tems as well as for comparison with CaFe2As2 in its 
collapsed tetragonal structure. We have calculated the 
electronic structure using 31x31x31 k-points in the full 
Brillouin zone and the resulting Fermi surface is shown 
in Fig.[2ja). During the calculation we relaxed the lattice 
position of the As atom, zas, from its experimental posi- 
tion of ZAs = 0.3665 to ZAs — 0.353(0) under the LSDA. 
The Fermi surface obtained is similar to that found in 
previous calculations for (Ba,Sr)Fe2As2^'^^. 

It is intriguing to note the similarity between this 
Fermi surface and that of MgB2 which also possesses two 
concentric warped cylinders at the zone corners as well as 
a flat 3D pocket in a similar position to the flared section 
of the dumbell surface. The primary difference is the lack 
of a further 3D pocket about F which occurs in MgB2, but 
not in CaFe2As2. In MgB2 it is the two outer cylinders 
which couple to the lattice to form superconductivityl^i. 



IV. THE 'COLLAPSED' TETRAGONAL PHASE 

In the collapsed tetragonal phase the crystal main- 
tains the lA/mmm symmetry of the high temperature 
tetragonal phase, but undergoes a dramatic reduction 
in the c-axis lattice parameter of approximately 6% 
while simultaneously the in-plane dimension increases 




FIG. 2: (Color online) a) The Fermi surface of CaFe2As2 in 
the high temperature tetragonal phase, b) The Fermi surface 
of CaFe2As2 in the 'collapsed' tetragonal phase at P=0.63 
GPcfcS. In contrast to the Fermi surface in the ambient pres- 
sure tetragonal phase we find just two Fermi surface sheets. 
One strongly 3D hole sheet that we refer to as the pillow as 
well as a strongly corrugated electron cylinder at the zone cor- 
ners. The 3D nature of the hole sheet makes the possibility 
of SDW instabilities due to nesting in the Fermi surface less 
clear. 



by approximately 2%. The unit cell is effectively 
squashed. At P=0.63 GPa the lattice parameters are 
a = 3.9780(1)A and b = 10.6073(7)Aa. We have cal- 
culated the electronic structure for this phase using a 
39x39x39 k-point grid. The resulting Fermi surface is 
shown in Fig. [DJb) . In Table U we also present the quan- 
tum oscillation frequencies and associated band masses 
for B parallel to c-axis. Experimental comparison with 
these values has the ability to verify the electronic struc- 
ture obtained in our calculations and the experimental 
mass enhancement can indicate the strength of correla- 
tions in the system. 

A critical observation of this investigation is the dra- 
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TABLE L Theoretical quantum oscillation frequencies and 
band masses for the collapsed tetragonal phase with the mag- 
netic field, B, parallel to the c-axis. 



Orbit" 


Frequency (kT) 


Band Mass irrie) 


a\ 


3.64 


1.25 


ai 


6.02 


2.17 


as 


7.19 


2.15 


/3 


16.0 


2.40 



"a denotes the warped cylinder around the zone corner in 
Fig. |2lb), /3 is the large pillow surface. 



niatic difference between ttie Fermi surfaces of the col- 
lapsed tetragonal phase and the high temperature tetrag- 
onal phase. The difference involves two key alterations. 
Firstly, in the collapsed phase we find a Fermi surface 
that is composed of just two sheets. Secondly, due to the 
further reduction in the c-axis lattice parameter the hole 
pocket is now very 3D. Therefore, the electronic struc- 
ture of this collapsed tetragonal phase is distinct from 
both the high temperature tetragonal phase and that of 
other members of the Fe-pnictide family. 

It has been suggested that magnetically mediated 
Cooper pairing is more robust in 2D systemsi^ and there- 
fore the 3D Fermi surface of the collapsed tetragonal 
phase may give rise to lower superconducting transi- 
tion temperatures than similar 2D systems. In fact if 
we consider the pressure induced superconducting tran- 
sitions of the (Ba,Sr,Eu)Fe2As2 family then we find sup- 
port for this idea. (Ba,Sr,Eu)Fe2As2 show maximal pres- 
sure induced transition temperatures of approximately 
32K, 28K and 29K respectively^'^. Each of these com- 
pounds shows a more 2D paramagnetic electronic struc- 
ture around the superconducting transition than we have 
found here for CaFe2As2. Therefore, the explanation for 
the lower transition temperature of 12K'^''', or absence of 
superconductivity as reported by Lee et al\ may rely on 
the different dimensionality of the Fermi surface in this 
collapsed tetragonal region. 

The dramatic change in Fermi surface topology in mov- 
ing through the transition from the high temperature to 
the collapsed tetragonal phase is produced by a reduc- 
tion in lattice volume of just 5%. This is indicative of 
a significant coupling between the electronic and lattice 
degrees of freedom. A phononic influence in combina- 
tion with magnetic interactions may therefore be active 
in producing superconductivity in these compounds as 
has been suggested for BaNi2As2^. 



V. THE ORTHORHOMBIC PHASE 

In the low temperature and pressure region of its 
phase diagram CaFe2As2 is antiferromagnetic (AFM) 
in the Fmmm orthorhombic structure. The lat- 
tice parameters are a=5.506(2)A, b=5.450(2)A, and 
c=11.664(6)A and zas = 0.36642(5)^8. The magnetic 



moment is 0.80{5)^b /Fe~^ . The conventional magnetic 
unit cell is the same as that of the crystal structure, but 
the symmetry is reduced due to the magnetic order. 

We have evaluated the electronic and magnetic struc- 
ture for the experimentally determined magnetic ordcr^^ 
using both the Generalized Gradient Approximation of 
Perdew-Burke-Ernzerhof (GGA) and the LSDA correla- 
tion functionals. For all calculations in this antiferro- 
magnetic phase we have used 23x23x10 k-points in the 
Brillouin zone. 

First, we applied the GGA to the AFM structure. Re- 
laxing the As position within the AFM calculation results 
in, ZAs = 0.3659, which matches the experimental value 
closely, but the predicted moment of 1.84/XB/Fe overesti- 
mates the experimental value by a factor of more than 2. 
Overestimates of this scale for the magnetic moment in 
Fe-pnictides while applying the GGA have been widely 
report edi^. 
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FIG. 3; (Color online) Plot of the force on the As atom(black) 
and the magnetic moment of Fe(red) versus the lattice pa- 
rameter ZAs in CaFe2As2. The forces and magnetic moments 
were determined under the Local Spin Density Approxima- 
tion. The main figure shows the region about Fz = where 
the As position is relaxed into its equilibrium structure. The 
inset is a view of this plot over a wider scale of zas- The 
solid horizontal line indicates where the force Fz ~ 0. The 
dashed lines in the main figure are to illustrate the value of 
the moment in /xs/Fe that occurs when Fz = 0. 

Secondly, we applied the LSDA to the antiferromag- 
netic structure while relaxing the value zas- We find 
that the relaxation of the As coordinate is significant, 
from ZAs — 0.3664 experimentally to zas = 0.3567. In 
Fig. [3] we have plotted both the force on the As atom, 
Fz, and the magnetic moment of Fe as a function of 
its relaxation to the plane of the Fe atoms. Very im- 
portantly, this figure shows that the moment approaches 
0.81 ± 0.02fiB/Fe as F^ tends towards zero ie. when the 
As is in its crystallographic equilibrium position. This is 
in close agreement with the experimental reported mo- 
ment of 0.80(5)^B/Fei^. To confirm this result we have 
applied the same method to SrFe2As2 and BaFe2As2. 
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TABLE II: Comparison between the magnetic moment of 
1,2,2 compounds predicted by the LSDA when the As position 
is relaxed in the antiferromagnetic phase and the experimen- 
tal valuesHi20|2i|22^ 



Compound 


LSDA 


ZAs 


Experimental 
{^iB)/Fe ZAs 


CaFe2As2 


0.81±0.02 


0.3567 


0.80(5) 0.36642(5) 
0.94(4)'' 0.3612(3) 


SrFe2As2 


0.97±0.03 


0.3507 


BaFe2As2 


0.86±0.02 


0.3444 


0.87(3) 0.35406(7) 



"The errors given are calculated from the convergence of the force 
on the As atom. 
''More recent measurementsSS have found /i = 1.01^s/-Fe. 




FIG. 4: (Color online) The Fermi surface of CaFe2As2 calcu- 
lated using the Local Spin Density Approximation with the 
As coordinate relaxed in the antiferromagnetic phase. 

The results are summarized in Table |lT] and give strong 
agreement with experiment. 

It has been suggested that overestimations of the mo- 
ment by Density Functional Theory may indicate quan- 
tum spin fluctuations. Such effects are not fully ac- 
counted for under Density Functional Theory which is 
mean-field in nature. Our study, however, suggests that 
a relaxation of the internal coordinates to the LSDA 
equilibrium position for As in the AFM phase is suffi- 
cient to permit the accurate prediction of the magnetic 
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moment. In some previous work the As coordinate has 
been relaxed in a non-magnetic calculation^!^ which has 
resulted in overestistimates of the moment of approxi- 
mately 1.5/is/Fe. Other work has shown the presence 
of a smaller moment under the LSDA in LaFeAsO^^ and 
BaFe2As2i^, but with less accurate correspondence with 
experiment than is shown in this study. Together, these 
findings suggest that the preceding overestimation of the 
moment should not be used as evidence for zero-point 
fiuctuations in this strongly antiferromagnetic material. 

The Fermi surface of CaFe2As2 for the LSDA cal- 
culation is shown in Fig. H) Most of the Fermi sur- 
face is gapped away by the additional ordering vector to 
leave the small pockets shown. Analysis of the associated 
quantum oscillation frequencies shows broad agreement 
with the recent experimental results of Harrison et al?^ 

The results presented here have the potential to mo- 
tivate further studies of the magnetic structure in these 
compounds that are both qualitatively and quantitatively 
accurate. 



VI. SUMMARY AND CONCLUSIONS 

The phase diagram of Fig. [1] illustrates that the 
CaFe2As2 Fermi surface changes dramatically through 
magnetic and structural phase transitions. This is sug- 
gestive that both magnetic and lattice interactions should 
be considered in the formulation of any model of super- 
conductivity in this compound. 

The greater three dimensionality of the collapsed phase 
Fermi surface suggest that pairing due to magnetic inter- 
actions may be weaker and this may explain the lower 
superconducting T,^, or absence of superconductivitjs^, 
in CaFe2As2 when compared to (Ba,Sr,Eu)Fe2As2 under 
pressure. 

Predictions for quantum oscillation experiments on the 
collapsed tetragonal phase have been given. We have also 
demonstrated that the magnetic moment in the AFM 
orthorhombic phase may be accurately calculated under 
the LSDA if the As coordinate is relaxed within an AFM 
calculation. 
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